APS/123-QED 



Carbon nanotube: a low-loss spin-current waveguide. 
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We demonstrate with a quantum-mechanical approach that carbon nanotubes are excellent spin- 
current waveguides and are able to carry information stored in a precessing magnetic moment for 
long distances with very little dispersion and with tunable degrees of attenuation. Pulsed magnetic 
excitations are predicted to travel with the nanotube Fermi velocity and are able to induce similar 
excitations in remote locations. Such an efficient way of transporting magnetic information suggests 
that nanotubes are promising candidates for memory devices with fast magnetization switchings. 

PACS numbers: 



The issue of magnetization dynamics in low- 
dimensional systems is currently one of the most studied 
in spintronics pj]. In particular, controlling how quickly 
the energy of a precessing magnetic moment propagates 
through wires is motivated by the possibility of build- 
ing miniaturized memory devices with fast magnetiza- 
tion switchings. This requires structures functioning as 
efficient spin-current waveguides, i.e., conduits capable 
of transporting spin information with little attenuation. 
Renowned for their high thermal and electronic conduc- 
tivities, for their long aspect ratios and coherence lengths, 
carbon nanotubes (NT) possess the ideal characteristics 
for acting as waveguides. NT have indeed been shown to 
function as conduits for electrons 0, Q and for phonons 
A natural question to ask, raised in this manuscript, 
is whether these materials are also efficient spin-current 
waveguides. 

Tserkovnyak Q proposed a mechanism for pump- 
ing spin current into a nonmagnetic metal caused by the 
precession of an adjacent magnetization. In this case, 
angular momentum from the moving magnetization is 
transferred to the conduction electrons, creating a spin 
disturbance that propagates throughout the metallic con- 
duit. This spin flow is produced without an applied volt- 
age, involves no net electrical current, and may be used 
to excite other magnetic units also in contact with the 
non-magnetic metal. Our strategy is to test this mecha- 
nism with NT in contact to local magnetic moments. 

Obstacles for identifying good spin waveguides lie in 
the difficulties of predicting the precise time evolution 
of the magnetization of a realistic system, taking into 
consideration the details of its electronic structure. This 
requires a full understanding of how a pulsed magnetic 
excitation evolves. Although previously treated semiclas- 
sically d, @ and quantum-mechanically 0, 13' ^ time- 
dependent quantum approach for this problem is still 
missing. In this manuscript we bridge this gap and pro- 
vide a quantum formalism that describes in time do- 
main how a localized magnetic perturbation propagates 
throughout a nanoscale system, more precisely a metal- 
lic NT. In doing so, we are able to show that NT are 
excellent spin-current waveguides. We demonstrate that 
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FIG. 1: (a) Spin susceptibility x™,o (in units of ft/eV) as a 
function of the excitation frequency oj (in THz) for a (4,4) NT. 
The distance between site m and the impurity is Xm = Soq. 
Inset shows NT with a single magnetic impurity. Arrow repre- 
sents where the electrons spin is probed, (b) Solid and dashed 
lines represent the x- and y-components of the spin distur- 
bance Sii{t); the thick solid line depicts [^^^(f)!. Spins are 
in units of hgnsho, where g is the electron spin g-factor, /is 
is the Bohr magneton and ho is the amplitude of the applied 
time-dependent transverse field; time is given in ps. Results 
were generated with a gaussian pulse centered at frequency 
uj = O.GTIfz and with a standard deviation a — 20.7 ps. 



they transport spin current across long distances in a 
dispersionless fashion and with little attenuation. More- 
over, we indicate that pulsed spin excitations propagate 
with a characteristic speed given by the electronic Fermi 
velocity and with a tunable degree of attenuation that 
depends on the excitation frequency. 

We consider a metallic NT doped with a single sub- 
stitutional magnetic impurity llo| , as depicted in the 
inset of figure 1. Assuming the existence of a localized 
magnetic moment pointing along the i-direction, we con- 
sider a perturbing time-dependent transverse magnetic 
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field h± (t) which sets the impurity magnetic moment into 
precession. Our treatment apphes to magnetic moments 
at temperatures above their Kondo temperatures which 
are nevertheless expected to be quite low in NT. We as- 
sume the Fermi-liquid representation and do not include 
any Luttinger-liquid effects. To determine whether a NT 
is a good conduit for spin currents, we must assess how 
such a localized magnetic excitation propagates across 
the structure. In other words, we must investigate how 
this excitation disturbs the spin balance of the system not 
only where the impurity is located but also, and more im- 
portantly, how the local spin dynamics is affected away 
from the impurity. This is naturally manifested by the 
spin susceptibility x, which reflects how the spin degrees 
of freedom of a system responds to a magnetic excitation. 

To calculate the spin susceptibility one needs the 
Hamiltonian describing the electronic structure of the 
unperturbed system, which we assume is given hy H ^ 

E<»j>,<t7»j ctcj<T + E<T(eono,<T + ^noafios) + Hz- 
Here, jij represents the electron hopping between nearest 

neighbor sites i and j, creates an electron with spin a 
in site i, eo is the atomic energy level of the magnetic im- 
purity occupying the site i = 0, noo- = cJ^cqo- is the cor- 
responding electronic occupation number operator, and 
U represents an effective on-site interaction between elec- 
trons in the magnetic site, which is neglected elsewhere. 
Finally, Hz plays the role of a local Zeeman interaction 



that defines the z-axis as the equilibrium direction of 
the magnetization. The Hamiltonian parameters can be 
obtained from density-functional-theory calculations so 
that the electronic structure of the doped system is well 
described [ll|, [l^] ■ Although we shall present results for a 
NT doped with Mn atoms, other substitutional magnetic 
impurities may be employed [lol |. In our calculations we 
fix the Fermi energy Ep = 0, and use jc,c = 2.7 eV, 
and ^Mn,c — 1-0 eV. We take the number of d-electrons 
in the Mn site no = l, eo = 0,C/ = 5 eV, and assume 
a local Zeeman energy splitting 5q = \ meV. Spin-orbit 
coupling is neglected for being very small compared to 
the electronic bandwidth of NT. 

The time-dependent transverse spin susceptibility is 
defined as XmjW = ~iQ{t){[S+{t). SJ {Q)]) . where 
Q{x) is the heaviside step function, and and S'~ are 
the spin raising and lowering operators at site m, respec- 
tively. The indices j and m refer to the locations where 
the field is applied and where the response is measured, 
respectively. In our case, a precession of the magnetic 
moment is induced at site j = 0, and we wish to observe 
the spin disturbance at an arbitrary site to. This re- 
sponse is fully described by Xmfl{t)- Within the random 
phase approximation, this susceptibility may be calcu- 
lated in frequency domain, and in matrix form it is given 
by x('^) ~ [^+X^{^) U]^^ where x*^ is the Hartree- 

Fock susceptibility, whose matrix elements are given by 



Here, j{io) represents the time Fourier transforms of 
the retarded single-particle propagators for an electron 
with spin a between sites m and j, and /(w) is the Fermi 
function. For a pristine NT 5^j(a;) may be analytically 
determined [l3| , and the integration in Eq. ([1} performed 
with great numerical accuracy. 

It is instructive to start by calculating Xm,o(w) given 
by Eq. (IlJ. A typical result is shown in Fig. la for a 
(4,4) armchair NT doped with a single Mn impurity. It 
depicts |Xm,o('^)| which is proportional to the amplitude 
of the spin disturbance at site ni due to a time-dependent 
transverse magnetic field applied at site 0. In Fig. la the 
site TO is a distance Xm = 3ao from the impurity, where 
oq « 2.46A represents the graphene lattice parameter. 
A very distinctive peak is evident in the w-dependent 
susceptibility, refiecting the existence of a resonance close 
to the Larmour frequency, as expected. 

Since we are interested in studying the propagation of 
magnetic signals, it is illustrative to look at the suscep- 
tibility as a function of time rather than frequency. The 
Fourier transform of the curve depicted in Fig. la gives 
the spin disturbance response to an instantaneous 5-like 



pulse h_L{t) applied at site 0. However, for practical rea- 
sons, we consider the response to a gaussian transverse 
magnetic field pulse given by h^{t) — Hq exp(— ap- 
plied along the x-direction (ho — hox). Here ho rep- 
resents the maximum field strength and a its standard 
deviation. Fig. lb shows how such a magnetic pulse ap- 
plied at the impurity location affects the spin balance of 
the NT at a nearby site {xm = 3ao). The thin solid and 
dashed lines represent the calculated spin components 
(S'i^(t)) and {S^{t)), respectively, and the modulating 
thick solid line is the magnitude of the transverse spin 
disturbance given by \{S^{t))\ = ^{S^,}^ + {S^^. To 
avoid congested figures the magnitude will be hereafter 
used to represent the spin dynamics on the NT. 

The fact that the electron spin probed on the NT dis- 
plays the precessional motion originated at the impurity 
is a clear indication of a flowing spin current emanating 
from the impurity into the NT. Earlier calculations have 
suggested that the existence of this current is sufficient to 
induce a dynamic magnetic coupling between dispersed 
impurities in NT, which is far more pronounced and long- 
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FIG. 2: Amplitude of the transverse spin disturbance (in units 
of hgnsho) as a function of time (in ps) at different distances 
Xm from the magnetic impurity. The curves labelled from (1) 
to (6) represent signals probed at different values of Xm ~ 
0.3,1.5,3,6,9, and 12 x lO^'ao, respectively. Inset plots how 
the corresponding probing positions Xm relate with the time 
r taken to reach a maximum spin disturbance. 



FIG. 3: Left axis depicts maximum amplitude of the pulsed 
excitation |{S'm)|max (in units of hgfiBho) probed at site m 
and is shown by the solid line with circular symbols. Right 
axis represents the width A (in ps) of the probed pulses shown 
by dashed lines with square symbols. Both quantities are plot- 
ted as a function of the separation Xm (in units of ao) from the 
impurity location. Panels (a) and (b) are for pulses centered 
at resonance and off-resonance frequencies, respectively. 



ranged than any other magnetic coupling of static nature 
0, S [H, [13) ■ To test how this precession propagates, 
in Figure 2 we show \{Sm{t))\ probed at different loca- 
tions on the NT. The six numbered lines correspond to 
spin disturbances probed at different sites along a line 
parallel to the NT axis containing the impurity site 0. 
The spin disturbance probed on other sites belonging to 
the same NT ring are virtually identical to the ones dis- 
played in Fig. 2. Small deviations are noticeable when 
the ring is very close to the impurity but disappear after 
a short propagation distance. This indicates that the ini- 
tial omni-directional excitation produced by the applied 
perturbation decays very fast into a cylindrical wave front 
that moves along the NT axial direction with a uniform 
speed. This propagation speed is found from the slope 
of the straight line in the inset of Fig. 2 showing the 
probing position Xm as a function of the time t spent by 
the pulse maximum to reach Xm- We find the propaga- 
tion velocity u = 1.4 x 10^ ni/s, which is precisely the NT 
Fermi velocity, indicating that the conduction electrons 
at Ep are the main carriers of the spin current. 

Figure 2 shows that, as the pulse moves, it is initially 
deformed but it quickly reaches a form that remains un- 
altered for asymptotically long distances. This can be 
seen when we focus on two features of the pulsed sig- 
nals, namely their maximum amplitude |(S'm)|max and 
their corresponding width, here represented by the mid- 
height separation A. Figure 3a plots |(S'm)|max and A as 
a function of the probing position Xm in the case of an 
excitation pulse whose frequency spectrum is centred at 
the same frequency as used for obtaining Fig. [T)d. Both 
quantities saturate after a short distance and are station- 
ary for Xm > 10* flQ. This is a remarkable result showing 
that the information contained in the magnetic pulse can 



be transported across very long distances without distor- 
tion. Surely, the asymptotic value for | (S'm) | max is not the 
same as the one probed near the impurity but it still has a 
sizable magnitude indicating that a considerable fraction 
of the energy contained in the magnetization precession 
can be used elsewhere. 

Furthermore, we can also control the fraction of the 
original pulse that reaches this asymptotic regime by 
selecting the central frequency in the excitation pulse. 
Figure 3b displays | (S'm)! max and A for the case of a 
pulse whose spectrum is centered at an off-resonance 
frequency cj — 0.24 THz. In this case the asymptotic 
value for | (Sm) |max is much reduced when compared with 
its counterpart in Fig. 3a. The existence of this non- 
decaying asymptotic contribution to the spin-current has 
been demonstrated recently for NT and atomic chains 
0, 13 but here we explicitly show that this can be actu- 
ally controlled by an appropriate selection of excitation 
frequencies. All the features presented so far indicate 
that NT can carry spin-current information contained 
within pulsed excitations without dispersing their form. 
This can be explained by the peculiar electronic struc- 
ture exhibited by NT, which is notoriously dispersionless 
at the Fermi level. In other words, extended states at the 
Fermi level have a linear dependence on the electronic 
wave vectors. In this case, group and phase velocities 
are interchangeable, meaning that every single frequency- 
component comprising a pulse will travel with the same 
speed, thus preventing any distortion of the pulse shape. 
If NT are to be used as spin-current waveguides, these 
features must be tested in the presence of structural dis- 
order. In fact, we have simulated the presence of isolated 
Boron substitutional impurities to explore their effect in 
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FIG. 4: Amplitude of the pulsed excitation as a function of 
time for the case of two magnetic impurities a long distance 
apart [xm = 6 x 10'*ao). Solid line shows the spin disturbance 
evaluated at the impurity where the excitation is induced 
whereas the dashed line is evaluated at the other impurity 
location. 



sured at the second impurity. Both lines have equally 
spaced peaks that are a time interval At — 2xm/v apart 
but are shifted by half that amount. This is easily ex- 
plained by the fact that precessing magnetic moments 
will always emit spin current to the surrounding conduc- 
tion electrons which will induce further precessions when 
interacting with other magnetic moments. Therefore, the 
induced precession probed at the origin appears as the 
first peak in Fig. |4l The subsequent pulse is seen to occur 
after the spin current has traveled all the way to excite 
the second impurity at site m. Once excited, this impu- 
rity precesses also emiting spin current, which will again 
induce another precession at the origin. This process is 
repeated indefinitely each time with smaller amplitudes. 
The amplitude reduction is inevitable because half of the 
energy stored in the pulsed precession will travel away 
from the other impurity. But the fact that the pulse 
produced at the origin reappears a long distance apart 
means that the energy stored in the precession of mag- 
netic moments can be noiselessly transported across the 
conduction electrons of NT, which confirms that these 
materials are ideal spin-current waveguides. 



the propagation of the spin disturbance pulse. We find 
that the pulse preserves its shape and amplitude, except 
very close to the B impurity, where it becomes slightly 
distorted. We notice that at larger distances, however, 
the pulse shape is restored as if no impurity is present. 
This suggests that the spin-current waveguide features 
reported so far are robust against the presence of non- 
magnetic disorder in NT. 

Finally, finding small values for the pulse amplitude at 
asymptotically large distances (see left axis of Figs. [3^ 
an[3]3) does not imply that the magnetic information is 
lost. To make this point more explicitly we add a second 
magnetic impurity a long distance from the original one 
[xm = 6 X lO^flo) and see whether the excitation pro- 
duced at the origin can be transported to the new loca- 
tion. Figure m shows the spin disturbance probed 
on both impurities. The solid line shows the magnetiza- 
tion precession at the origin whereas the dashed line rep- 
resents the spin disturbance as a function of time mea- 



In summary, we have shown that metallic NT are ex- 
cellent spin-current waveguides and are able to trans- 
port magnetic information across long distances with 
minimum dispersion and with very little loss. Spin dis- 
turbances induced by localized magnetic excitations are 
shown to propagate throughout the length of a metal- 
lic NT with speeds of the order of 10^ m/s. Moreover, 
we show that the energy stored in a magnetic preces- 
sion can be used elsewhere when the spin current trav- 
eling through NT interacts with other magnetic objects. 
These features turn NT into ideal components for fast- 
response memory devices. The experimental testing of 
our predictions requires the generation and detection of 
spin currents. One possibility is to excite a spin with 
a laser pulse, and monitor a second spin located some- 
where downstream with an STM tip. Since the STM 
current can be modulated by a coherent spin precession, 
one should be able to observe the second spin being ex- 
cited by the spin current in the NT. 
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